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Abstract 

We present the result for the finite part of the two-loop sunrise integral with unequal masses 
in four space-time dimensions in terms of the C)(s*’)-part and the 0(s^)-part of the sunrise 
integral around two space-time dimensions. The latter two integrals are given in terms of 
elliptic generalisations of Clausen and Glaisher functions. Interesting aspects of the result for 
the 0(s')-part of the sunrise integral around two space-time dimensions are the occurrence 
of depth two elliptic objects and the weights of the individual terms. 


1 Introduction 


The two-loop sunrise integral with non-vanishing internal masses is the first and simplest inte¬ 
gral in quantum field theory, whieh cannot be expressed in terms of multiple poly logarithms. It 
has already received considerable attention in the literature [HHHl . However, the question which 
generalisations of multiple polylogarithms appear in the evaluation of the two-loop sunrise inte¬ 
gral is still open. A first and partial answer has been given for the two-loop sunrise integral in 
two space-time dimensions [|l3lll5l . In two space-time dimensions the two-loop sunrise integral 
is finite and does not require regularisation. Furthermore only one graph polynomial contributes 
in the Feynman parameter representation. It turns out that in two space-time dimensions the two- 
loop sunrise integral is given as a product of a period with a sum of elliptic dilogarithms. This has 
been shown for the equal mass case in lfT3l and for the general unequal mass case in llTSl . The 
arguments of the elliptic dilogarithms have a nice geometric interpretation as the images in the 
Jacobi uniformisation of the intersection points of the variety defined by the graph polynomial 
with the integration region. 

In this paper we study the two-loop sunrise integral around four space-time dimensions. Our 
motivations are as follows: First of all, four space-time dimensions is the physical value and 
it is the sunrise integral in four space-time dimensions which enters precision calculations in 
high-energy particle physics. Secondly, we hope to learn more about elliptic generalisations of 
poly logarithms. In particular we are interested in generalisations with higher weight or depth. 
Various elliptic generalisations of (multiple) polylogarithms have been discussed in the mathe¬ 
matical literature [fT3l[T5l[T9] - l2^ . The definitions of the generalised functions differ among the 
articles in the mathematical literature and depend on the postulated mathematical properties one 
would like to keep in going from the (multiple) polylogarithms towards the elliptic case. We are 
interested to find out, which definition nature has chosen in quantum field theory. 

Unlike in two space-time dimensions, the sunrise integral in four space-time dimensions is 
divergent and requires regularisation. Dimensional regularisation is the method of choice. In 
D = A — space-time dimensions, the £-expansion of the two-loop sunrise integral starts at 
1 /e^, corresponding to two ultraviolet divergences. The terms proportional to 1 /e^ and 1 /e are 
very well known and our main interest in this article are the terms proportional to E°. 

Using dimensional-shift relations we can relate the 0(E°)-term around four space-time di¬ 
mensions to the 0(E®)-term and the 0(E^)-term of the sunrise integral around two space-time 
dimensions. In the equal mass case the coefficient of the 0(E^)-term in this relation vanishes 
and the 0(£®)-term around four space-time dimensions is related to the 0(£®)-term around two 
space-time dimensions alone. In this paper we treat the more general case of unequal masses and 
discuss as a specialisation the case of equal masses. 

The 0(£®)-term of the sunrise integral around two space-time dimensions is already known 
[|13[I15II . In this paper we compute the 0(E^)-term of the sunrise integral around two space-time 
dimensions and express it in terms of generalisations of Clausen and Glaisher functions towards 
the elliptic case. We briefly recall that the Clausen and Glaisher functions are related to the real 
and imaginary parts of the classical poly logarithms. Generalisations of multiple poly logarithms 
towards the elliptic case start to make their appearance in physics ffT3l[T5l[24ll25ll . The 0(£^)- 
term of the sunrise integral around two space-time dimensions gives us information on elliptic 
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generalisations of multiple polylogarithms of depth greater than one. In addition we observe that 
this term is not of uniform weight, it eontains terms of weight three and four. We diseuss in detail 
the oeeurrenee of the weight four terms. 

This paper is organised as follows: In seetion|2]we introduee the two-loop sunrise integral. 
Seetion[3]is devoted to variables assoeiated to this integral. Our prineipal method to ealeulate the 
integral is based on differential equations, and the differential equations obeyed by the two-loop 
sunrise integral are diseussed in seetion HI The solution of any differential equation requires 
boundary values, and the boundary values for the two-loop sunrise integral are presented in 
seetion [51 In seetion we introduee generalisations of the Clausen and Glaisher funetions to 
the elliptie setting. These funetions will be needed to express our final result for the two-loop 
sunrise integral. In section|7]we give the result for the 0(£^)-part of the sunrise integral around 
two space-time dimensions. Using dimensional-shift identities we can relate the terms up to the 
finite part in the £-expansion of the two-loop sunrise integral around four space-time dimensions 
to the 0 (£®)-part and the 0 (£^)-part of the sunrise integral around two space-time dimensions. 
This is done in section [ 8 l Section |9] discusses the two-loop sunrise integral around four space- 
time dimensions for the special case, where all masses are equal. Section [T0| is devoted to a 
detailed discussion of the transcendental weights occurring in our result for the 0 (£^)-part of the 
sunrise integral around two space-time dimensions. Finally, section [TT] contains our conclusions. 
In an appendix we included explicit formulae, which are too long to be presented in the main 
text of this paper. 


2 Definition of the sunrise integral 


The two-loop integral corresponding to the sunrise graph with arbitrary masses is given in D- 
dimensional Minkowski space by 


W1V2V3 




nil,m2, 


y) = 


( 1 ) 




v-D f d^k\ d^k 2 




D 

in 2 


D 

in 2 


k^ + m 


2)V: 


-k? -|- m 




{p-k\ -k2)^ + m2 


V 3 • 


In eq. ([U) the three internal masses are denoted by mi, m 2 and m 3 . Without loss of generality we 
assume that the masses are ordered as 


0 < mi < m 2 < m 3 . (2) 

The arbitrary scale p is introduced to keep the integral dimensionless. We denote by v = Vi -f 
V 2 -I- V 3 the sum of the exponents of the propagators. The quantity denotes the momentum 
squared (with respect to the Minkowski metric) and we will write 

t = (3) 

We perform our calculation in the region defined by 

-r > 0, (4) 
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and in the vicinity of the equal mass point. Our result can be continued analytically to all other 
regions of interest by Feynman’s iO-prescription, where we substitute t —)■ t + iO. The symbol 
+i0 denotes an infinitesimal positive imaginary part. 

Where it is not essential we will suppress the dependence on the masses m,- and the scale 
/j and simply write 5viV2V3(f^T) instead of In terms of Feynman 

parameters the two-loop integral is given by 


W 1 V 2 V 3 


(D,t) = 


r(v-D) 


r(vi)r(v 2 )r(v 3 ) 






VI —1 V 2 —1 V 3 —1 






r 


D 


(0 


(5) 


with the two Feynman graph polynomials 


U = XlX2+X2X'i+X'iX\, 

‘J = —xiX2XT,t-\-{xim\-\-X2m2-\-XT,m^)U. ( 6 ) 

The differential two-form (0 is given by 

(0 = xidx2f\dx2,—x2dxif\dx2,+x2,dx\f\dx2- (7) 

The integration is over 

a = {[xi : X2 : jcs] G > 0 ,/= 1,2,3} . ( 8 ) 


In order to facilitate a comparison with results in the literature we remark that the definition of 
the sunrise integral in eq. Q is in Minkowski space. In a space with Euclidean signature one 
defines the two-loop sunrise integral as 

^ Tl^ + {K^ + miy^(^{P-Ki-K2f + miy' 

The momenta in Euclidean space are denoted by capital letters, while the ones in Minkowski 
space are denoted by lower case letters. We have with P^ ~ —p^ the relation 

= ‘^ViV2V3 • ( 10 ) 

We are primarily concerned with the integral Si 11 (4 — 2£, t). This integral has a Laurent expan¬ 
sion in £, starting at e 


Siii(4-2£,t) = e 


-2ye 


■^s[ j j ^ (4, t) -1- -s[ j ^ ^ (4, t) -b s[ y {4,t) + 0 (£) 


( 11 ) 


_n'\ _ 

The first two terms Sj j j (4, t) and S\ j(4, t) are rather simple, while the term Sj"^} (4, t) is in the 
focus of this paper. With the help of dimensional-shift relations we relate the integral Sm (4 — 
2 £, t) to the integral Si 11 (2 — 2£, t), which has an expansion in £ starting at £®: 




2 ( 0 ) 


Siii( 2-2£,0 = si‘!}(2,0+esSO(2,t) + C>(£2) 


:( 1 ) 


( 12 ) 
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The term 5|°j(2,r) has been given in ifTSll in terms of elliptie dilogarithms, the term 
will be given in this paper. 

In using dimensional-shift relations we will encounter simpler integrals, obtained by pinching 
internal propagators. For the sunrise integral any integral obtained by pinching propagators is a 
product of tadpole integrals. The tadpole integral is given by 


Tv (D,m2,^ 


D 

"T 


d^k 


D 

in 2 


[—k^ + ni^) 


r(v-f) 

r(v) 


m 




(13) 


3 Variables related to the sunrise integral 

The sunrise integral defined in eq. O depends for a given space-time dimension D on the vari¬ 
ables t, m\, ni^, m 3 and It is clear from the definition that the sunrise integral will not change 
the value under a simultaneous rescaling of all five quantities. This implies that the integral de¬ 
pends only on the four dimensionless ratios t / m| / and m^h^. It will be convenient 

to view the non-trivial part of the integral as a function of five new variables 

222 

mfm^m, 

q, wi, W 2 , W 3 and -(14) 

The variables wi, W 2 and W 3 satisfy 


W 1 W 2 W 3 = 1, (15) 

therefore there are again only four independent variables. The variables q, wi, W 2 and W 3 are 
closely related to the elliptic curve defined by fF = 0 . 

Before we discuss this change of variables we first introduce some convenient abbreviations. 
Let us start by denoting the pseudo-thresholds by 


/ri = mi-(-m 2 — m 3 , /r 2 = mi — m 2 -|-m 3 , /r 3 = —mi-|-m 2 -|-m 3 , (16) 


and the threshold by 


/i4 = mi-|-m2-|-m3. (17) 

In the vicinity of the equal mass point we can assume m 3 < m\ + m 2 , or equivalently /ri > 0. 
With this assumption and with the ordering of eq. Q we have 

0 < A'l < A^2 < A'3 < A'4- (18) 


We further set 




(19) 
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The variable D will appear in the algebraie prefactor of our result. It will be convenient to 
introduce the monomial symmetric polynomials in the variables m^, and These 

are defined by 


= I (20) 

o 

where the sum is over all distinct permutations of A few examples are 


Mioo 

= m\ + ni 2 + m 2 , 


Mni 

222 
= m^m2m3. 


M210 

42 42 42 42 42 42 

= mjm2+ m2m3+m3mj+m2mj+m3m2+ mim3. 

(21) 


In addition, we introduce the abbreviations 

ooo ooo ooo 

A = /Ji/J2)U3/r4, 5i = — mj+m 2 + m 3 , 62 = mj—+ 63 = m^+m 2 —m 3 (22) 

and 


vi = 


1 -e/v^ 

^ ' 5i 


I — i 




V2 = 


i + '-f 


1 + / 


;VA 


I — i 




V3 = 


I — i 


'-ST "“‘■ST 

In the vicinity of the equal mass point we have 

A>0, 5, >0, /€ {1,2,3}. 




(23) 




(24) 


The variables vi, V 2 and V 3 are then complex numbers of unit norm. In addition they satisfy 


V 1 V 2 V 3 = 1. (25) 

We now describe in detail the change of variables from the set t, m{, m^, mg,^^ to the set specified 
in eq. (fT4l) . The equation 


IT = 0 (26) 

is a polynomial equation of degree three in the Feynman parameters xi, X 2 and X 3 . The cubic 
curve intersects the domain of integration in the three points 

Fi = [l:0:0], F 2 = [0 : 1 : 0], F 3 = [0:0:l]. (27) 

The cubic curve together with a choice of a point O G {/’i,/’ 2 ,^ 3 } as origin defines an elliptic 
curve. All three choices will lead to the same Weierstrass normal form, given by 

= 4{x — ei){x — e 2 ){x — e 2 ), with ei+e 2 + ^ 3 = 0 , (28) 
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and 


= 24^4 + 2Mioo? + A + , 

^2 = (^-?^ + 2Mioo?+ A-3\/Dj , 

^3 = ~ 4Mioo? - 2A) . 


The modulus k and the eomplementary modulus k' of the elliptie eurve are given by 


k=J^^, k' = 
V ^1-^2 


The periods of the elliptie eurve ean be taken as 

.Jdx 


' gi -g3 
gi — g 2 


¥1=2/— = (fc), ¥2 = 2/— = k'), 

2 y Di jym 


y D? 


y D? 


^3 ^3 

8^/^ f{x-e 2 )dx 4^2 g 4 i^(^x-e 2 )dx 4z>2 

(^i = — / - = —{K{k)-E{k)), ^2 = — -= ^£(k). 

J y J y 

ei ei 

K{x) and E (;c) denote the eomplete elliptie integral of the first kind and seeond kind, respeetively: 




= / dti 


l-x^t^ 
l-t 2 ■ 


The Legendre relation for the periods reads 


¥i4'2-¥24'1 = -^• 


The Wronskian is given by 


d d 4-4 (3^2-2Mioot + A) 

We denote the ratio of the two periods \\f 2 and \|/i by 


and the nome by 


q = e 


1 



The nome q will be one of our new variables, and we can think of the variable q as replacing the 
variable t. A useful relation for expressing q in terms of t (or vice versa) is given by 


_^^_ 

(a^i - 0 - 0 - 0 {^4 - 0 n (x)^* 

In this equation, Dede ki nd’s rj-function, defined by 


ri(x) 


oo oo 

n (1 =<?^ n 

n=l n=l 



(37) 


(38) 


appears. The left-hand side of eq. (l37l) has a power series in t, the right-hand side of eq. (1371) has 
a power series in q. Eq. (1371) can be used to express t as a power series in q, or vice versa. We 
have 


t 


A 2 


+ 0{q^)^ 


q = —t 




A 2 


■ + 0 (t 2 ), 


(39) 


therefore t ^ 0 implies (? —)■ 0 and vice versa. It remains to define the variables wi, W 2 and W 3 . 
These are given by 



P, = 7t 


F{ui,k) 

K{k) 


Ui = 


gl — g2 
^j,k ~ g2 


Xj,k = e3 + 


m]ml 
/ ■ 


(40) 


In the definition of ut we used the convention that {i,j,k) is a permutation of (1,2,3). In the 
definition of P, the incomplete elliptic integral of the first kind appears, defined by 



dt 

y/(l -t2) (I 


In the case <? = 0 (or equivalently t = 0) we have 


(41) 


lim w i 

q-^O 


Vj, ye {1,2,3}, 


(42) 


while in the equal mass case mi = m 2 = m 3 we have 


271/ 


Wi 


J \mi=m2=mi 


= e^, ye {1,2,3} 


(43) 


for all q (or equivalently all t). There is a simple geometric interpretation for the variables wi, 
W 2 and W 3 : We recall that the points Pi, P 2 and P 3 are the intersection points of = 0 with the 
integration region a. One of these points is chosen as origin of the elliptic curve. The set 


|wi,W2,W3,Wi 


,W 2 ^ 



(44) 


is obtained from the images in the Jacobi uniformisation of the two points not chosen as origin 
by considering all three choices of origins. More details can be found in ifTSl . 
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4 Differential equations 


In this section we present the (fourth-order) differential equation for the sunrise integral 5iii(D,t) 
in D space-time dimensions with arbitrary masses. From this equation we deduce a differen¬ 
tial equation for the 0(£^)-piece The latter differential equation is then solved up to 

quadrature. 

In two space-time dimensions the integral Sm (2, t) = (2, t) satisfies a second-order dif¬ 

ferential equation ffTTI : 






(45) 


The coefficients p 2 , p\ and po as well as ps are collected in appendix lAl The left hand side 
defines a Picard-Fuchs operator 


4°^ (2) = P2^ +Pi^^+P0- 


d 

'dt^ dt 

The periods and \\f 2 are solutions of the homogeneous differential equation IfT^ 


(46) 


= 0, /G{1,2}. 


(47) 


We will encounter several differential operators. We will use the notation 

7^4(7)), (48) 

where r denotes the order of the differential operator and D denotes the associated space-time 
dimension. In the case where D is an integer and a superscript j is present, this superscript 
denotes the order in the £-expansion to which this operator belongs. Finally, i is a label to 
distinguish differential operators with identical r, D, j. 

In D dimensions the integral 5iii(D,t) satisfies a fourth-order differential equation. 


d^ d^ d^ d 
+7^3 "7^ +7^-7^ +7^1^ +7h 
dr dr dr dt 


>5111 (7),0 


P^[ci2Tl2+C\2,T\2,+C2'iT2'i\- (49) 


Here we used the abbreviation 


T,j = (50) 

The explicit expressions for the coefficients P 4 , P 3 , P 2 , Pi, Po and C12, C13, C23 are rather long 
and given in appendix El There are several possibilities to obtain the fourth-order differential 
equation: It can either be derived by using the relations given in flbl, by using the program 
“Reduze” [I26[I27II or by the algorithm given in [(281. Eq. (l49l) defines a fourth-order Picard-Fuchs 
operator 

d^ d^ d^ d 

+ + (51) 
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The Picard-Fuchs operator L^{D) has a polynomial dependenee on the number of spaee-time 
dimensions D. Around D = 2 — 2£ we can write 


U{2-2t) = 

j=o 


(52) 


Of particular relevance for the 0(£^)-part (2, t) will be the operators (2) and (2). The 

operator (2) factorises 

if (2) = 4"] (2) 4">(2) if (2). (53) 

The differential operator L^\2) is the one we already encountered in eq. (|4^ . The two other 
factors t[^](2 ) and t[^](2 ) are first-order differential operators: 




^ 2 (St^-\-2MiQQt-\-ld 

^ (15^2 - 2MiQQt - 3A) (3^2 - 2Mioot + A) ^ 

60^2 — 12A/ioot^ ~ ( 6 OA /200 ~ 88 A/ 110 ) t — \2 M\qqA. 
(15^2 - 2Mioot - 3A) (3^2 - 2 Mmt + A)^ 

(15^2 - 2Mioot - 3A) ^ - (30t - 2Mioo) ■ 


Solutions to the homogeneous equations 


4 °] ( 2 ) ¥«(0 = 0 , Lfl{2)Mf,{t)=Q 


(54) 


(55) 


are 

Vfl (0 ~ (3t^ — 2Mioot + a) -|-2A/ioot + 7A) , 

Mfb{t) = 15t2-2Mioot-3A. (56) 

Substituting the £-expansion of 5iii(2 — 2£,t) given in eq. (fT^ into the D-dimensional differ¬ 
ential equation (l49l) gives a coupled system of differential equations for 5[23j(2,t), where the 

differential equation for 5jy(2,t) will involve the lower order integrals 5y\(2,t) with i < j. 
This system can be solved order by order in £. At order £® one finds 

4°i(2) Lfl{2) if {2) sS°l(2,t) = -32A/V(l5t2+14Mioot + 77A). (57) 

From eq. (1451) we know already that 

Lf (2) 5S‘;i(2,0 = (58) 

and eq. (1571) reduces to 

Lfl{2) Lfl{2) p^{t) = -32^V(l5t2+14Mioot + 77A), (59) 
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which is easily verified. 


At order we have 


L<">( 2 ) 4 ">( 2 ) X< 1 ',( 2 ,<) = /,(<), 


(60) 


with 

/i(<) = 


- 4 ‘'( 2 ) X<“',( 2 ,<) 

{912^4 + 1344Mioot^ + (9088Miio-4416M2oo)2‘^ + 512MiooAt+ 112A2 


m? . T T TN. m? 


+d m\, m3) In —A + J (4 m\, m3, m\) \n—j+d m\, m^) In 


a /2 


(61) 


and 


d (t, mj, m 2 , m 3 ) = —320?"^ — (352mf + 212m^ + 212m^ t’ 

+ (l440m| + 1744m2 + 1744m3 — 3120mjm2 — 3120mjm3 — 36\6m2m^^ 

+ (^-5AAm\ + 212 ml + 212 ml + 1296mV2 + 1296mVI - 1024mfm^ - 1024mfm^ 

—212m\m^ — 272m2m3) t 


+ (22Am\ — 112m2 — 112m3 — 112mjm2 — 112mjm3 + 224m2m3) A. 


(62) 


Eq. (t60l) is a fourth-order differential equation for the integral j (2, t). The fourth-order differ¬ 
ential operator in eq. (l60l) faetorises into two first-order differential operators and a seeond-order 
differential operator. The inhomogeneous term h {t) has a Taylor expansion in t, starting with 
Eq. (1^ is easily solved for L2*^^(2) (2,t): 


( 0 ) 


Lr( 2 ) X< 1 ',( 2 ,<) = / 2 (<), 


with 


h [t) = 


(63) 


(64) 


i 

Ci^ifbit)+C 2 ^b{t) J 


'^aih)dti 


+^bit) 


\\faitl)dti 


n 

I 


h {h)dt 2 


7 Pl,bitl)\\fb{tl) Pl,ait 2 )^a{t 2 )' 


where C\ and C 2 are two integration eonstants to be determined from the boundary eonditions. 
In eq. (1631) we used the notation that 


Lfl{ 2 )=pi^bj^+pQ^b- 


( 65 ) 


Eor the homogeneous solutions we have 

C\'^b{t) = Cl (l5t^ — 2A/ioot — 3A) , 
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t 

Ci^bit) [ . = ^(3?2 + 6Mioo?-7A)l (66) 

0 

The differential equation (l63l) is now of the same type as eq. (1^ . only the inhomogeneous term 
differs. It ean be solved with the same methods as used for eq. (1451) by ehanging variables from 
t to the nome q. This ehange of variables is described in detail in lITSl . One finds 


MuP'O = 


C3V[ti+C4V[t2- 


n 


dqi 

<?l 




dqi h{q 2 )^\{q 2 ) 

<?2 np2{q2)W [qiY 


(67) 


C 3 and C 4 are two further integration constants, which need to be determined from boundary 
conditions. Eq. (1671) gives the solution for the integral 5[|j(2,t) up to quadrature. We still need 

to determine the integration constants. In addition we would like to express the integral 5^ j j (2, t) 
in terms of elliptic generalisations of multiple poly logarithms. Let us also mention that eq. (1671) 
can be used to obtain the (^r-expansion (or equivalently the t-expansion) of sj j j ( 2 , t) to high orders 
in q (or t). 


5 Boundary values 


In this section we give the boundary values at t = 0 for the first two coefficients of the £- 
expansion of 5iii(2 — 2 £,t) (i.e. 5[*jj(2,0) and 5[jj(2,0)) and for the first three coefficients 
of the £-expansion of 5iii(4 — 2£,t) (i.e. 5[^j^^(4,0), 5[^j^^(4,0) and (4,0)) We will use the 

boundary values 5j|j(2,0) and Sj°^(4,0) together with regularity conditions at t = 0 to fix the 
integration constants C\, C 2 , C 3 and C 4 of the previous section. 

Let us start with the boundary value at t = 0 of the sunrise integral in D = 2 — 2£ dimensions. 
We have 


5iii(2-2£,0) 


r(l + 2e)(r)‘+"‘y' 

a 


(D 

(^xim\ + X2rn2 + x^rn^ 


( 68 ) 


By a change of variables we can relate this integral to the one-loop three-point function in 4 -f 2£ 
space-time dimensions (please note the sign of the £-part) with massless internal lines and three 
external masses. The change of variables can be found in lfT 2 l and the result of the one-loop three 
point function can be taken from 02911301 . One obtains for the £-expansion of Sm (2 — 2£,0) 


'£eJs[{\{2,0) = e^^r{l + 2e) 

j=o 


j - \ 1 2£ ry 


P 


1 r(i+£)- , . , 

2£2r(l+2£) £ 


,(69) 


with 


// = 


-Vj) ^ 27^1 (-2£,-£;1-£;v;) - (-V- 2 ^^! (-2£,-£; 1-£;v ■ ^ 


.(70) 
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The expansion of the hypergeometrie funetion reads 

2Fi(-2£,-£;1-£;x) = l + le^Lhix)+e^2Lh{x)-4U2,i{x,l)] +0{e^). (71) 

We obtain for (2,0) and j ^ (2,0) 


2 ^- 


2 3 


*^111(250) — ^^2/ Li 2 

.2 3 




*^111 (2,0) = % E ^ {- 2 Li 2 ,i {vj, 1) -Lig (v;) + 2 Li 2 ,i (v^- \ l) +U 3 (v ■ ^ 


\/A “i2i 


j 

-2 In 



Li 2 (v;) -Li 2 


(72) 


The definition of the multiple polylogarithms is given in eq. (1791) . Let us now turn to the values 
of the eoeffieients of the Laurent expansion of Siii(4 — 2£,t) at t = 0. At t = 0 we find 

Mioo 


.(- 2 ) 

*111 


(4,0) = 


2 ’ 


( 11 " (4.0) 


2 H 

3Mioo , rn\ 

H—^In 


' 2 \ 2 
mf \ mi 

+^ln 
. A' / A' 


n (4,0) = (2,0)- + 3 



The expression for (4,0) is obtained from the expression given in |[ 6 l by noting that the 
function 


L = Li 2 
1 

+- 


m 2 

-/ 

mi 


mi 


A3 + Li 2 -A3 + ^2 + T A3 


m 2 


1 


In A 3 + 


m 2 

mi 


- In A 3 + 


mi 

m 2 


2 

3, m 
T-ln 


mi 


m\ 

In^, 

mi 


A 3 = 


1 


m3 — mf — m2 + V — 


2 m\m2 

can be written in a more symmetrical way as 
1 


^ 2 L 


Li2 (vi)+Li2(v2)+Li2(v3)-Li2 (vj M -Li2 (V 2 M -Li2 (vg ^ 


(74) 


(75) 
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The constants of integration Ci, C 2 , C 3 and C 4 are determined as follows: In section[ 8 ]we discuss 
dimensional-shift relations. These relations allow us to relate the integral 5iii(4 — 2£,t) to the 
integral 5iii(2 — 2£,t), and we can determine the integration constants from the following four 
conditions: 

1. The requirement that 5^ j (4, t) is regular at t = 0, i.e. there is no pole at t = 0 for Sj*] ^ (4, t). 

2. The boundary value 5^ j (4,0). 

3. The requirement that j j (2, t) is regular at t = 0, i.e. there is no logarithmic singularity at 

t = Oforsi5^i(2,t). 

4. The boundary value 5^ j j (2,0). 

The explicit expressions for C\ and C 2 are given in appendix El Condition 3 implies C 4 = 0. With 

¥i(? = 0) = (76) 

one finds for C 3 

C 3 = 2^M‘il(2.0). (77) 

6 Generalisations of the Clausen and Glaisher functions 

In this section we introduce elliptic generalisations of the Clausen and Glaisher functions. These 
generalisations will show up in our final result. 

The classical polylogarithms are defined by 

00 j 

= L (78) 

,/=i 7 

and the multiple polylogarithms by 

°° il-l fyt-l-l Jl J 2 Jk 

Li„i,„2,...,„,(xi,X2,...,Xfc) = £ ^ (79) 

,/i = 172=1 7A:=1 7 l 72 Jk 

The sum representation gives rise to a quasi-shuffle product and one has for example 

Lijil ('l'l)Lin 2 (' 1 - 2 ) ~ Lini,n 2 (■^l;-^ 2 ) “l“Lin2,7Ji ■hLi„j ^„2 (' 1-1 ' 2 C 2 ') ■ (80) 

We recall that the Clausen functions are given by 

( ^ [Li„ -Li„ , n even, 

Cl„((p) = < (81) 

[ 5 [Li„ -fLi„ , n odd, 
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while the Glaisher functions are given by 


( 2 [Li„ + Li„ (e "P)] neven, 

Gl„((p) = I (82) 

[ Yi [Li„ (e'‘P) - Li„ , n odd. 

The Clausen and Glaisher functions correspond to the real and imaginary part of Li„(e'^), and 
the actual assignment depends on whether n is even or odd. 

Let us now consider the elliptic setting. In IfTSl we considered the following generalisation 
depending on three variables jc, y, q and two (integer) indices n, m: 

oo oo j k 

ELi„;„(x;y;^) = (^3) 

We define the weight of ELi„;m(v:;y;( 2 ') to be w = n + m. The definition is symmetric under the 
exchange of the pair (x, n) with (y, m). The two summations are coupled through the variable 
q. In the special case q = \ the two summations decouple and we obtain a product of classical 
poly logarithms: 


EEi;j;OT (.^?y? 1) — Lin (x) Ei,jj (y). 


(84) 


In addition we introduce the following linear combinations 


En;m(^5y5^) (85) 

{ i [^Un{x)-^U„{x^^) +ELi„;„(v:;y;^)-ELin;m(x^^;y^;<?)], n + meven, 

jhin (x) + jLin ^) + ELi„;„ {x; y,q)+ ELin;„ (jc^ ^; y^ ^, n + m odd. 

The special case {n,m) = (2,0) appeared already in IfTSl . Eq. (I85l) gives the generalisation to 
arbitrary indices {n,m). In general, the functions En-m{yy',q) are not symmetric under the ex¬ 
change of the pair {x, n) with (y, m), nor do they have for m 7 ^ 0 a uniform weight. The functions 
En-,m{yy\q) can be thought of as elliptic generalisations of the Clausen and Glaisher functions. 
In particular we have 


limEi;o (e'‘‘’;y;<?) 

= Cli((p), 


limE 2;0 (e'‘‘’;y;<?) 

= Cl2((p), 


limEs;! (e'‘P;y;<2') 

= Gl3((p). 

(86) 


We now turn to the multi-variable case. In order to keep the notation simple in this case, we 
introduce a prefactor c„ and a sign 5 „, both depending on an index n by 


c» = i|(i+0 + (i-0(-i)"] 


I, n even, 
z, n odd. 


sn^{-\r = 


1 , n even, 
— I, zz odd. 


(87) 
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With the help of these definitions we ean write the definition of the funetions uni¬ 

formly as 

^n-m{x\y\q) = ( 88 ) 

^n+m ^ 2 ^^” ^ ELin;;?; (x ,y ,q^ 

It turns out, that in addition to the funetions defined above we only need a single depth two 
elliptie objeet in order to express all results of this paper. This funetion depends on five variables 
^ 1 , ^ 2 , yi, yi, q and five (integer) indiees n^, n 2 , m\, m 2 and o. This new funetion is defined as 
follows: 



Cn+m 


2 Lin ('*f) T ELin;m (^x,y,q') 


^ni,n2',mi,m2',2o{2Ci,X2,yi,y2^q') — 


(89) 


Crii+mi Cfi 2 -\-m 2 


ELifjj(xj ,yi, ^) 5n[+miELi,2j-|-o;mi+o 


Lin2 {-^ 2 } ^n2+m2^^n2 i^2 


^2 L 


CX3 CX3 CX3 CX3 

+ E E E E 

jl = lfcl = l72 = 1^2 = l 


jl kt —j\ —ki 

^lEl -^ni+mi -^1 El 


Jl k 2 ~jl —k 2 

X^y^ - Sn 2 +m 2^2 y 2 


q 


jlki+j 2 k 2 


n\ 


JTk' 


ni 1/112 




( 71 ^ 1 + 72 ^ 2 )" 


Please note that this definition is asymmetrie in the quadruplets (Ri,mi,xi,yi) and (n 2 ,«^ 2 , 2 C 2 ,y 2 )- 
Let us briefly diseuss the weights of the individual pieees. The first term, eonsisting of produets 
of the form 


ELi,2j_^2>;mi+o ('’f'l jEI 5 Lin2 (-^ 2 ) 


(90) 


is of weight wi = ri -f- ^2 + mi -|- 2 o, while the quadruple sum of the form 


CX3 oo oo oo 


E E E E 

71 = 1*1 = 172 = 1^2 = 1 




(91) 


is of weight W 2 = Ri + «2 + mi -|- m 2 -|- 2o. For m 2 = 0 the two weights eoineide. 

For o > 0 we ean express this funetion as an o-fold iterated integral over q (the variables xi, 
X 2 , yi and y 2 are treated as eonstants in the integration): 


Enj ,n2;mi,m2:2o (2fl 5 X 2 ,yi ,y25 — 



dq 2 

q 2 



(•Ifl jEI 5 <7o) ^n\\mi 


(•ii:i;yi; 0 )] 

^n 2 \m 2 {.x 2 -,y 2 ',qo) ■ 


(92) 


Ineq. (I92l) the asymmetry withrespeet to the quadruplets (ni,mi,ji;i,yi) and (n 2 ,m 2 ,X 2 ,>' 2 ) mani¬ 
fests itself by the faet, that the eonstant part with respeet to the Taylor expansion in q is subtraeted 
out from E„j;mj(xi;yi;( 5 i), but not from E„ 2 ;m 2 (-i^ 2 ;E 2 ;<?)- This subtraetion ensures that the inte¬ 
grand has a Taylor expansion starting at < 3 '^ Therefore the integral is well defined at the lower 
integration boundary ^ = 0 and no regularisation is needed. 
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Let us write an explicit example. At depth two there is only one combination of indices, 
which will appear in the results of this paper. The indices are given by (ni,n 2 ) = (0,1), = 

(—2,0) and 2o = 4. We have 


Eo, 1;-2,0;4 {xi,X2\yi,y2\-q) = 

^'^ELi2;o(.vi;yi ;-<?)-ELi2;0 ““S') ^ \ Eil (jC2)+Eii (.^2 ^ 


OO OO oo CXD 

+ E E E E 


*? 


21 = 1 /ti = l 22=U2=1 72(71^1 +72^2) 

The corresponding integral representation reads 


-’fi'jl’ j {^2y^ +-^2“^^E2*^ j (^_qyih+i2h 


(93) 


hH 1 


/ 


Eo, 1;-2,0;4 (-^1 A2;E1 ,E2; = 



0 0 


Eo ;-2 (-"Cl ;ei ; 0)] Ei;o (.r2;E2; -qi) ■ 


(94) 


7 The result for (2, f) 

In this section we present the 0(£^)-term (2, t) of the sunrise integral around two space-time 

dimensions in terms of elliptic generalisations of the Clausen and Glaisher functions discussed 
in the previous section. 

Eet us start by recalling that the 0(£®)-term 5[*Jj(2,t) of the £-expansion of 5iii(2 — 2£,t) is 
given by 


4°! ( 2 .<) 



(95) 


with 


£(0) = E2 ;o(wi;- 1;-<?)-I-E2 ;o(w 2;-1;-<?)+E2 ;o(v £3 • (96) 


This motivates the following ansatz for y (2, t): 




7t 


(97) 


We find that is given by 

e(i) = 


2 ^ 
^7=1 


2 ^ 

rrij 


6Ei,o(-l;l;-^) + 

7=1 L 


Ei;0 { wf ,\',- q ) { wf ,-\-,- q ) 


,£( 0 ) 
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Y- {Li2,i (wj, 1) -Li2,i (^Wj ^ l) +Li3 (wj)-Li3 (wj ^ 


+3 In (2) 


Li2 {wj) - U 2 



3 

+ Y, [ 4 Eo,i;_ 2 ,o; 4 (>Vy,W;;-l,-l;-^) - 6 Eo,i;_ 2 , 0;4 {wj,Wj;l,-\;-q) 

7=1 

+6Eo,i;-2,0;4 
3 3 

+ L L [2Eo.i;-- 2,0;4 -Eo,i;-2,0;4(vvji,Wy2;-l,-l;-<2') 

71=172=1 

-Eo, 1 ;- 2 , 0;4 -q )] 

3 

+ J^E 3;1 {wj-,-\-,-q). 

7=1 


(98) 


Eq. (I98l) is one of the main results of this paper. A few remarks are in order: Eq. (I98l) gives the 
result for entirely in the variables q,wi,W 2 ,W 3 and m\m^m^l. Although not as elegant as 
the result for in eq. (|9^ . it is still far from trivial that ean be expressed in a few lines. 
As arguments of the (elliptie) multiple poly logarithms only the variables and the values ±1 
occur. 

The first line gives all terms proportional to the 0(£°) result These terms are all of 
weight one, multiplying which is of weight two, yielding the total weight three. The sec¬ 
ond and third line contain terms, which depend on the variables Wj, but not explicitly on q. These 

are ordinary multiple poly logarithms. Terms like these are expected, since 5^ j ^ (2, t) depends not 
only on the graph polynomial IT, but also on the graph polynomial U. The ordinary multiple 
polylogarithms are all of weight three. The next four lines contain functions of depth two. Sur¬ 
prisingly, we only encounter the function Eo,i;- 2 , 0;47 which is of weight 1 — 2 -(- 4 = 3. The last 
line gives probably the most interesting part: It contains the elliptic poly logarithm E 3 ;i, an ellip¬ 
tic generalisation of the Glaisher function GI 3 according to eq. (I 8 ^ . The elliptic polylogarithm 
E 3 ;i is not homogeneous in the weight, having parts of weight 3 and parts of weight 4. The 
occurrence of E 3 ;i is discussed in more detail in section [TO 

Eurthermore, the result shows an explicit ln(2). These ln(2)-terms are spurious and cancel 
with ln( 2 )-terms from Ei;o(— 1 ; 1 ; —^) and EQj-^2,o-A{wj-i “ 1 ; “ 1 , 


8 The result for the sunrise integral around four space-time 
dimensions 

In this section we discuss dimensional-shift relations, which relate the integrals 5iii(4 —2£,7) 
and 5i 11 (2 — 2£, 7 ). This allows us to express the finite part sj*] j (4, t) of the sunrise integral around 
four space-time dimensions in terms of the already calculated integrals Syj(2,7) and 5[jj(2,7). 
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The divergent parts and are rather simple and listed for eompleteness. 

There are several relations, whieh relate integrals in D spaee-time dimensions to integrals in 
(D + 2) spaee-time dimensions. The Tarasov relations read II31[I32II 


5'viV2V3 (^T) ~ ViV2*S'(vj_|_1)(v 2-(-1)V3 + 2, t) -|-V2V35vj (v2 + l)(V3-t 1) T 

+ViV35(vj_^i)v2(v3+i) (T> + 2,t). (99) 

From eq. dS]) we find 

A^^^5'viV2V3 (£>T) = ViV2V35(vj + 1 )(v 2 + 1 )(V 3 + 1 ) (^ + 2,t) ■ (100) 

For a fixed spaee-time dimension we ean always reduce any integral to a linear combination of 
sunrise master integrals and simpler integrals. Inverting the relations, which relate integrals in 
D dimensions to integrals in (D 4 - 2) dimensions, we can express Si 11 (4 — 2£, t) in terms of four 
sunrise master integrals for D = 2 — 2e and simpler tadpole integrals. A possible basis of sunrise 
master integrals is given by 

Siii(D,t), ;/2 ^Siii(D,0, /^Siii(D,0, (101) 

In this basis one finds for Sy ^ (4, t) 

^m(4T) = (2) sS;^i(2,t) + lLf (2) sffi(2,t)+7?)°), (102) 

where the remainder function contains the contributions from the tadpoles. The differential 

operator L 3 ^^( 2 ) factorises again: 

(2) = l[-^\2) Lf (2), (103) 

with l[ ( 2 ) given by 


—33t^ -\- 13A/ioot^ + (22A/110 — 15A/200 ) ^ ~ 3A/iooA d 
2At (3t^ — 2M\QQt A) (5t^ -|- 2M\QQt -1- 7A) dt 
69t^ — 22A/ioot+ 9 A /200 ~ 104^110 
^ 24t (3^2 - 2MiQQt + A) (5^2 + 2Mioot + 7A)' 


The explicit expressions for the differential operator (2) and the remainder function i?)®) are 
rather long and not provided here. We have used eq. (11021) internally for the determination of the 
integration constants C\ and C 2 . 

In order to present our final result on sj*] ^ (4, t) it is better to use an alternative basis of sunrise 
master integrals given by 


5iii(Z),t), (AO- 


dm\ 


dm^ 


dm] 


(105) 
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This is the basis of sunrise master integrals used in [O. We recall that 








{D,t) 

= S222{D + 2,t), 

' ^ 2*^111 
omj 

= —S2ii{D,t), 

' ^ 2*^111 
am^ 

= —Sui{D,t), 

' ^ 2*^111 
am^ 

= —Sn2{D,t). 


Let us define for G { — 1,0} two partial differential operators 


(106) 


if (2) = cf + £c; 

/=1 


W 2 

m 


* dm} 


(107) 


For k — —\ the coefficients are given by 


C 


C'- 


(- 1 ) 




[(2m^ — mj — m|) t — 3mj — 3m^ — mjmj — m^m| + , 


- 1 ) 


= i(Ct 




(108) 


while for k = 0 we have 

[—Mioot^ + (—26M2oo + 36Miio)t+M3oo + 8M2io —78Miii] 


= 


[-t^ + (I2mj - 3m) - 3ml) t • 
+66mym^] . 


m,- 


16m{ — 16mf — 3m} m} — 3m} m} 


(109) 


Please note that all coefficients of (2) vanish in the equal mass case. For j (4, t) we find 


(4.,) = if(2.,) +if (2)5ffi (2.,) + f 


( 110 ) 


- E 

{iJ,k)eZj, 


4mf + m) + ml) t + mf Imj + ml) — 2m)ml 




m? — 2m) — 2m^ t d-m} (m) + m^ 


2m)ml 



6/i^r 

2t^ — (24ml d- 6tn) + 6m|') t + 2m) — m^ — m) — 6m? (m) + m)^ + 12m?m| 


I2^^t 


■In 
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The notation (/, 7 , k) 6 Z 3 stands for a sum over the three eyelie permutations of (1,2,3). Eq. (IllOl) 
expresses the finite part j (4, t) of the sunrise integral in 4 — 2£ dimensions in terms of Sj j ^ (2, t) 
and sf[\ (2, t), derivatives thereof and simpler terms. The results for j ^ (2, t) and ^ (2, t) have 
been given in section |71 Please note that L 3 (2) vanishes in the equal mass case. Therefore it 
follows that in the equal mass case the 0(£)-part sj (2, t) does not affect (4, t). However, in 

the unequal mass case the 0(£)-part 5[jj(2,t) is required for 5|°^(4,t). Eq. (IllOl) together with 
eq. (|9^ is the main result of this paper. 

Eor completeness we list also the divergent terms of 5i 11 (4 — 2£, t). These read 


5n?(4,<) 

■^ill 


Mm 

t 3Mioo fmj 


( 111 ) 


9 The equal mass case 


This section is devoted to the equal mass case mi = m 2 = m 3 = m. In the equal mass case the 
result for the sunrise integral around four space-time dimensions is significantly simpler. We 
have already seen that the contribution from 5[jj(2,t) in eq. (IllOl) drops out in the equal mass 
case. Eurthermore we have now a second-order differential equation for all D Q: 

lltiyt — 9m^) {t — m^) -|- [3 (4 — D) -|- 10 (D — 6) tm^ -|- 9Dm"^] — 

-f (D - 3) [(D - 4) t + (D + 4) m2] } Si 11 (D, t) = -3{D- if [Ti {D)f, (112) 

with Ti{D) = Ti{D,m^,/j^). A convenient basis of sunrise master integrals is therefore 


Siii(D,t), f/—SniiD,t). 

In terms of this basis the dimensional-shift relation simplifies and we have 


5'iii (4 —2£,t) = 


1 


(t-(-3m2] (t —m2] (t —9m2) d 


6(1-2£)(1-3£)(2-3£) 

+ 


{t — {t — 9m2] {t'^ llrn^t — 87m^) 

1“ £ A 






1 ^ dt 

5'iii (2 —2£,t) 


(113) 


‘S'lii (2 —2£,t) 


-|-3 (1 — £) 


67,2 + 

T ^ A 


m^ 


nr 


[ri(4-2£)]2 . 


Eor the coefficients of the £-expansion we have now 
3m2 


5'llf(4,<) = 


2 n 


2 ’ 


(114) 


(115) 
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9m^ 3m^ 

f^ln 


4/j2 2 ^ /j 


( 


m 


‘^lll(4’0 — 0..2 


13? 9m^ 3m^ 




C 2 -U -2 


8/J^ jA^ 2fA 

(t-m^) {t-9n?) ( 0 ) 




9nP' \ , / 

In I ^ 


3m^ ^1 






cl-; p I (^ + 3 m 2 ) (?-m 2 ) (?- 9 m 2 ) J (q) 

12/ l2/ 

again showing the absence of 5^ j j (2, t) in the expression for 5^ j (4, ?). 


10 Weights 


In this section we discuss in more detail the transcendental weights associated with individual 
terms and in particular the occurrence of the functions £ 3 ;!. 

The building blocks of all our expressions are the multiple polylogarithms of the form 


°° fl^l 7/1-1“ 1 y.71 yil yjk 

Li„i,„2,...,n,(jCl,JC2,-,-^lt) = L L ••• L 

71 = 172=1 7it=l ^2 Jk 


the generalisation of the classical polylogarithms to the elliptic case 


00 00 ^j\ k\ 

ELi„j;„j (xi;yi;/ = L L 

n=\ki=i h H 


and quadruple sums of the form 

00 00 00 00 Jl/l ^jlki+jiki 

V V V V ^1 -^1 -^'2 yi ^ _ 

hhhh {hh+hhr 


(116) 


(117) 


(118) 


For xi = X 2 — ...Xk — 1, yi =yi = 1 and i? = 1 the summand is in all cases a homogeneous 
function of the summation variables. We defined the transcendental weight as the negative of the 
degree of homogeneity with respect to the summation variables in this case. Thus eq. (11161) is of 
weight wi = «! + ^2 + ■■■ + kik, eq. (11171) is of weight W 2 = n\-\-m\, and eq. (11181) is of weight 

W3 = Hi + ^2 + mi + m2 + 2o. 

With this weight counting, the function 


ELi3;i (xi;yi;< 2 ') (119) 

is of weight 4 and we would like to discuss the occurrence of this function in the result for 
s| j ^ (2, t). It turns out, that the occurrence of this function can be related to the fact that 

In^ (120) 
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should be counted as weight 2. The latter fact can be seen as follows: We have for example 

y^ELi„;„ (/;>';<?) Jinx' = ELi„+i;„ , (121) 

0 


but 


q 

y^EEi„;„ (x;>';(2'') Jin*?' = mAn+i-m+i{x\y,q). (122) 

0 


Thus an integration with respect io d\nq = dq/q increases the weight by 2. The final result for 
( 2 , t) does not contain any In^-terms, however a convenient method for calculating ( 2 , t) 
splits the integral into two parts, each part containing In^-terms. The In^-terms cancel in the 
sum, but leave the functions EEi 3 ;i as remainders. 

In order to understand how this happens, it is sufficient to consider the equal mass case 
mi = m 2 = m 3 = m for sj (2, t). In the equal mass case the differential equations for j j (2, t) 
(and ( 2 , t)) simplify and we find with 


.( 0 ) 


d^ 

2 ,equal ^ ^ 


d 


/^2,equal = t [t - n?) [t - 9n?) , Pl,equal 


(123) 

3t^ —20tm^ + 9m^, po,equal = t —3m^, 


the following differential equations 

r(0) , 

^2,equal"- 


p(®) 5^®) (2 t) — 


-6p2, 


(124) 


(2T) = + 




-3t^ + lOtm^ + 9m^)-3t + 5m^ 

' dt 


do) 

’ill 


( 2 , 0 - 


In order to find 5[jj(2,t), let us make the ansatz that 5[jj(2,t) consists of a part proportional to 


( 2 , t) and a remainder 5^ j j ( 2 , t): 


^ln(2,») = di’i(2,<)-l-fiWdn(2,»). 


(125) 


The differential equation for 5[j ^ (2, t) is then 


4')l,„d'l)(2,l) = 12Pln(^^)+6PF,(,)-[(^: 


d'^Fi{t) 


2p2,equal - lOtm^ - 9m'^^ ^ 


.V / dF\{t) 2 

Tp2,equal 1“ Pi,equal 1 \-3t 5m 


dt 


M 

’ll! 


( 2 , 0 - 


(126) 


We are free to choose Fiit) (different choices for F\{t) will lead to different expressions for 
5[jj(2,t)) and a convenient choice for Fi{t) will be the one, which eliminates dsf^^{2d)/dt on 
the right-hand side of eq. (11261) . This amounts to the choice 

1 1 1 

Yt~ 


d —3P'+\0tm^ + 9m^ 

dt ^ 2t(t — m?-) (t — 9trP-) 


trP- t 


9m?-' 


(127) 
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and hence 


m(t — m^) — 9m^) 

3/^ 


fi(<) = 



(128) 


It turns out that this choice does not only eliminate the terms proportional to the derivative of 
but also the terms proportional to itself and the differential equation for 

5[y(2,t) reads 

^lequal*^!!! equal (0 i f2,equal (0 = ^lu ' (^29) 

The inhomogeneous term of this differential equation is significantly simpler than the generic 
inhomogeneous term of eq. (11261) . But please note that Fi{t) and /i,equal (0 contain both terms 
proportional to ln(t) (which can be related to ln(—( 5 r)-terms). 

In order to determine 5[jj(2,t) from eq. (11291) one needs first the homogeneous solutions. 
These are spanned by tjti and \\f 2 , defined in eq. (l3TI) . For the case at hand it will be convenient 
to use instead of the basis {t|/i, V(t 2 } the basis given by 

\\fu v^iln(-^). (130) 


We then write the full solution as 


with S 


Sn\i2j) = ciV|ti+C2Viln(-^)+5i5\^^pg^i^i(2T) 

111 special ^ F^om the boundary values one finds now 

3 


(131) 


Cl = 


C2 = 


Ini 


: iLi3(r3)-Li3 (rg M - 2 Li2i (r3,1)+Li3 (r3) -Li2i (rg ^ 1) - Li3 (rg ^ 


In(3) (^Li2(r3)-Li2 (^ 3 ^))}’ 
^ ^ '^Li2(r3)-Li2 


2 27tz 

In the following we will denote by 


(132) 


{2m\ 

r, = exp(^—J 


(133) 


the p-th root of unity. For special (^’0 Hi^ds along the lines leading to eq. (l67l) 


.special 


:( 1 ) 


00 00 


111.special _ 9; 

j=lk=\ 




r dqi r dq 2 

J q\ J qi 

0 0 


^ 3 -^ 3 ^ 


fc2(-l)^(-^2)^'^/2,equal(-?2). (134) 


24 














In order to perform the integration we need the ^-expansion of / 2 ,equal- The ^-expansion of / 2 ,equal 
is given by 


^2,equal (<?) — 31n ( (J') 


(135) 


CX3 OO -j 

+ 12 EE 7 

j=\k=\ J 


6 r ^2 + ^3 + ^ 3 '^ ■ 


6r^- 


6 r. 


?j 


+ 3r^ + 3r^T 


Irir, — IrVr. — Ir'ti — 2r 


7; 


12 


12 


12 


.11^' 

12 


q- 


jk 


Of partieular relevanee for the weight 4 part is the term 31n(—^). We have 


2^1 ^ f dqi f dci2 

n 2ip^^ J qi J q2 

; ik-Iq q 


(r^'-r3^)fc2(-l)^(-^2)^'^ln(-^2) = 


(^ELi2;o(r3;-l;-^) -ELi2;o (^3 j 

(^EEi3;i (r3;-l;-^)-EEi3;i j . 


(136) 


Counting the weight of ln(—^) as two, we notiee that in eq. (11361) all terms are of weight 4. The 
two terms on the right-hand side of eq. (11361) eombine nieely with parts of the boundary terms to 
give 


-^^ln(-<5r)E2;o(r3;-l;-<?) + 3^E3;i(r3;-l;-^). (137) 

2 7t 7t 

We note that eq. (11371) eontains terms of weight 3 and 4. The term proportional to ln(—(?) in 
eq. (11371) eaneels the eorresponding logarithmie singularity of ln(—( 2 -) (or equivalently ln(t)) in 
Fi(t)5yj(2,t). The seeond term proportional E 3 ;i explains the oeeurrenee of E 3 ;i in the final 
result for 5^ j j (2, t) in eq. (1981) . 


11 Conclusions 

In this paper we presented the result for the (7(£^)-part of the sunrise integral around two spaee- 
time dimensions. The result is expressed in terms of generalisations of the Clausen and Glaisher 
funetions towards the elliptie ease. The (7(£^)-part gives us information on elliptie generalisa¬ 
tions of multiple polylogarithms of depth greater than one. It is worth noting that the (7(£^)-part 
of the sunrise integral around two spaee-time dimensions eontains terms of weight three and four. 
It is not of uniform weight. We diseussed in detail the oeeurrenee of the weight four terms. Using 
dimensional-shift relations we expressed the finite part of the sunrise integral around four spaee- 
time dimensions in terms of the 0(£°)-part and the (7(£^)-part of the sunrise integral around two 
spaee-time dimensions. 
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A The coefficients of the differential equation 


We first recall from ifTTlI the coefficients of the second-order differential equation for Siii(2,t). 
This differential equation reads 


(f- 


d 




5iii(2,t) = 


(138) 


where po, pi P 2 and pj, are given by 

P2 = t{t-pfj{t-pl){t-nl){t-pl){3t^-2Mioot + ^), 

Pi = 9t^-32Mioot^ + (37M2oo + 70Miio)t4-(8M3oo + 56M2io + 144Miii)t3 (139) 

— ( 13 A /400 — 36 M 310 -I-46A/220 ~ I 24 A/ 211 ) 

— (— Sfffsoo T 2437410 — 1637320 — 9637311 -1- 1443722i) t 

— (47600 ~ 63/510 -I- 153/420 — 203/330 -I- 1837411 — 123/321 — 63/222), 

Po = 3t^—737ioot'^ +( 23 / 200 +1637iio)t^-I-( 63/300 — 143 / 210 ) 

— (537400 — 837310 + 637220 — 8372ii) t + (37500 — 3374io + 237320 + 837311 — 1037221) • 


The polynomial p^ appearing in the inhomogeneous part is given by 


P3 


-2{3t^-2Mioot + Af (140) 

2 2 2 

fn^ tYt‘~ fft 

+2c(t,mi,m2,m3)In—2 -f 2c(t,m2,m3,mi)In^ + 2c(t,m3,mi,m2)In^, 


with 


c(t,mi,m2,m3) = 

(—2mf + m2 + m3) + (6m| — 3m2 — 3m3 — Irr^tr^ — lm\m^ + 14m2m3) 

+ 6mi + 3m2 + 3m3 + 1 \m\m\ + 1 \m\m^ — %ni\m\ — Sm^m^ — 3m\ni^ — 3ni^m\ 

+ ^2mi — — 5 m\tr^ — 5 m\m^ + m^m^ + m^m^ + Am\n^ + 

-\-3m\m\ + 3m\m^ — bm^m^ + 2m\m\m^ — rn\m\m^ — m\n^m^ . 

In D dimensions the integral 5iii(D,t) satisfies a fourth-order differential equation: 


^ d^ ^ d^ ^ d 


( 141 ) 


^lll {D,t) = p^ [ci2Ti2 + Ci3ri3 + C23723] . (142) 


The coefficients Pj (with 0<j< 4)read 

P 4 = 8^3 {t - pI) {t - pI) [t - pI) {t -p^ 4 ) [{1 - D) - 2 {D- 3) 37ioot + (13 - 3D) A] , 
P 3 = 4^21 (SD^ - 7 ID + 252) - {2D^ - 82D + 500) Mmt^ 
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- (33^2 - 343D + 520) M 2 oo?^ - (14Z)2 - 130D - 144) Mnot^ 

+ {52D^ - lllD + 2040) - {lOD^ - lOOD + 312) M 2 io?^ 

- (152^2 - 2552D + 9744) Mi11(13^^ - 503D +1780) 

+ (36^2 - 684D + 2192) + (124^2 - 1332D + 4592) M2nt^ 

- (46^2 - 362D + 824) M22ot^ - (1 8D2 + 46D - 508) M^QQt 

+ (54^2 + 138D - 1524) M^wt - (120^2 + 904D - 6288) M^nt 
+ ( 132d2 + 1532D - 9528) M22it - (36^2 + 92D - 1016) M32o? 

-3D(3D-13)a2}, 

^2 = 2 ? I - (9D2 - 183^2 + 1224D - 2688) - {%D^ - 62D^ - 326D + 2360) Mioot^ 

+ (5 - 999^2 + 5718D - 9360) M 2 oo?^ - {6D^ -1SD^ + 828D - 3696) Mi 

- (24^2 - 852^2 + 6948D - 15360) Mgoo?^ + (56^2 - S52D^ + 4516D - 8928) Mzio?^ 
+ (192^2 - 4392^2 + 31560D - 71856) Mi 11 ?^ - (19^^ - 13D^ - 1636D + 6640) 

+ (12^2 + 28^2 - 2304D + 8528) + (20^^ + 112D^ - 9536D + 27248) M2ii?2 

+ (14^2 - 202^2 + 1336D - 3776) M 22 o?^ - (l 14^2 - 414D - 312) M^oQt 
+ (342^2 - 1242D-936) M 4 io/‘+ (144^2 -2472^2+ 7608D+ 1872) Ms 1 + 

- (288^2 - 4260^2 + 12732D + 1872) M22\t - {22SD^ - 828D - 624) M32ot 
-3D{D-2) (3D-13)a2}, 

= (VD^ - 175^2 + 1610^2 - 6440D + 9408) 

+ (14^4 - 270^2 + 1708^2 - 3792D + 1440) Mioot^ 

- (27^4 - 767^2 + 7182^2 - 27352D + 36480) M2oot^ 

+ (38^4 - 782^2 + 6172^2 - 22384D + 31296) Mnot^ 

- {12D^ + 52D^ - 2712^2 + 16688D - 29760) M^qqP 

- (20^4 - 532D^ + 4888^2 - 19376D + 28608) M2wt^ 

- (56^4 - 2168^2 + 24688^2 - 112576D + 179136) Mnit^ 

+ (17^4 - 273^2 + 1342^2 - 1416D - 2880) M^oot^ 

- (20^4 - 180^2 + 184^2 + 2592D - 7488) 

- (UOD^ - 2220D^ + 11656^2 - 19680D - 1728) M2nt^ 

+ (6^4 + 186^2 - 2316^2 + 8016D - 9216) M22ot^ 

- {2D‘^ - 34^2 + 292^2 - 1088D + 1248) M^qqI 
+ (6D4 - 102^2 + 876^2 - 3264D + 3744) M^i^t 

+ (104^4- 1128^2 + 2992^2+ 1440D-7488) M 31 + 

- (220^4 - 2460^2 + 7736^2 - 3648D - 7488) M22it 


27 



- (4D4 - 6SD^ + 584^2 - 2176D + 2496) 4 / 320 ? 

-D {D -2){D- 4) {3D - 13) A^, 

Po = (D - 3) (D - 4) I - (D^ - 2+ 146D - 336) - {3D^ - 5ID^ + 258D - 360) Mioo?^ 

+ {2D^ - lOD^ + 56SD - 1280) M 200 P - {\2D^ - 204-D^ + 1144D - 2192) Mi 10 ?^ 

+ {6D^ - 66D^ + \20D + 240) Msoo?^ - {6D^ - 50D^ + 40D + 336) M2io?^ 

- (28^3 - 300D^ + 584D + 1392) {D^ - 33D^ + 182D - 240) M 400 ? 

+ (12^3 - 148^2 + 560D - 624) M 310 ? - (12^^ - 84^^ + 48D + 144) M 211 ? 

- {22D^ - 230D^ + 156D - 768) M 220 ? 

+ (D - 2) (3D - 13) [(D -4)M3oo -{D- 4)M2io + (2D - 16)Miii] A}. (143) 

In the inhomogeneous term we have 
C 12 = (D - 2f (D - 4) { (D^ - 13D + 42) 

+D^ (2mj + 2m2 + 4m3) P — D {24m\ + 24m2 + 40m3) + {62m\ + 62m2 + 76m3) 

+D^ (—4mi — 4m2 + 2m3 + \0m\m\ + lOm^m^ + Sinfm^) 

+D (62mi + 62m2 + 22m3 — \32m\m^ — \32m\m^ — 108mim2) 

+ (-198mt - 198m^ - 180mj + 450m?m^ + 450m^m^ + 348m?m^) 

+D^ 2mj — 2m\ — Am\ + 2m\ni^ + 2ni\m\ — ^m\m^ — %m\ni^ + 14mim3 + XAni^m^ 

+ 16mjm2m3) t 

+32Dm3 {3m\ + 3m2 + 2m3 — 5m\m^ — 5m^m^ — 4mim2) t 

+ (^2m\ + 42m2 — 180m3 — 42mjm2 — 42mjm2 — 2%%m\rr^ — 2SSm\rn^ + 426mjm3 
+426m2m3 + 192mfm2m3) t 

+ (3D — 13) [D (—m| — + m 3 + 2m^m2) + 4mj + 4m2 — 6 m 3 + 2trP^n^ + 2n^m^ 
- 8 m?m^] A}. (144) 

ci 3 and C 23 are obtained by permutation of the masses. 

B The integration constants 

The explieit expressions for the two integration eonstants Ci and C 2 appearing in eq. (1631) are 
given by 


Cl 


4 

3 


(M300 - M210 + 5 Mi 1 1 ) Sj®) (2,0) + ^2 2 ^ 
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- (imj — mj — — m^mj — In I ^ 


— mj — — Sm^mj + Im^m^ 


— — Sni^ml — 


— m3 — mfm2 — mjm3 + Irr^tr^ In^ 

— m| — m2m3 — m^mf + Irr^rr^ In^ ^ ^ 





{m\+m\- 

~ {m\^m\- 


lm\ + mfm3 + m^m^ — 


lm\ + m^mj + m^mj — 


2 m 2 + m 3 m 2 + mjm^ — Im^ni^ 





, (145) 


C2 = (7M600 - 44M510 + 113M420 - 77M411 + I2IM321 - 152M330 - 258M222) (2,0) 


+/<^ ^(121M3 oo-121 M2 io + 270Miii) 


+ ( 2 mj 




■ ml) In^ f 2^ j + (2m3 —ml — m^) In^ ( 


ml\ , f ml\ T T .ox. f ml\, f mV 


+2 [ml + ml- 2ml) In j In [j^j + 2 (m^ + - 2ml) In In 


'mV\ . f mV 


+2 [ml + ml- 2ml) In j In j 
2 r 

32m}®+ 55 [ml + ml)m\— (355ml + SSSm^-400mlml) m^ 

+5 (ml + ml) (SSml + SSm^ — 124m2m3) m| — (ml — ml)^ (l45m2 + 145m3 + 546mlml) m\ 

—2(ml — ml)‘^(ml + ml) In 
2 r 

+ 32m2® + 55 (m 3 +mj) m^ — (355m3 + 355mj — 400m3mi) m® 

+5 (m 3 +ml) (83m3 + 83mj — I24mlml) m\ — (ml — ml)^ (I 45 m 3 + 145m| + 546m3mj) m^ 


-2(ml-miy 


Dhn n 
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7A2 L' 


32mf + 55 (mf + mg ml - {355m\ + 355m^ - 400mfmg 


mo 


+5 (mj +mg (83mj + 83m2 — 124mjmg m 3 — {m\ — mg^ (l45mi + 145m2 + 546mjmg 


mo 


-2K-mg"K + mg 


In 


(n^ 

W 


(146) 
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